Multimode quantum interference of photons in multiport integrated devices 
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We report the first demonstration of quantum interference in multimode interference (MMI) 
devices and a new complete characterization technique that can be applied to any photonic device 
that removes the need for phase stable measurements. MMI devices provide a compact and robust 
realization of NxM optical circuits, which will dramatically reduce the complexity and increase the 
functionality of future generations of quantum photonic circuits. 
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Photonics is a leading approach to realizing future 
quantum technologies that promise secure communica- 
tion pQ, tremendous computational power [2] and the ul- 
timate precision measurements [3]. Optical waveguide 
circuits on silicon chips have demonstrated high lev- 
els of miniaturisation and performance @H7], however, 
the construction of large-scale multi-port devices from 
2x2 directional couplers makes them unwieldy and sen- 
sitive to wavelength. Multimode interference (MMI) de- 
vices promise a straightforward implementation of robust 
multi-port circuits, however, multi-mode operation could 
prevent or perturb quantum interference, which to date 
has not been demonstrated in this architecture. Here 
we show that quantum interference can be realised in 
MMI devices, observing a quantum interference visibil- 
ity of V = 95.6 ± 0.9% in a 2 x 2 MMI coupler. We 
further demonstrate operation of a 4 x 4 port MMI de- 
vice with photon pairs, which exhibits complex quantum 
interference behaviour. We have developed a new tech- 
nique to fully characterise such multi-port devices based 
on measuring quantum interference for all possible two 
photon input and output combinations. This approach 
removes the need for phase sensitive measurements and 
may find applications for a wide range of photonic de- 
vices. These results show that MMI devices can operate 
in the quantum regime with high fidelity and promise 
substantial simplification and concatenation of photonic 
quantum circuits. 

Quantum technologies aim to harness superposition 
and entanglement to enhance communication security, 
provide exponential computational advantage for partic- 
ular tasks [8] , including factoring [9] , database search [10] 
and simulation of important quantum systems [TT] , and 
reach the ultimate limits of precision in measurement [3 
and lithography [T2] . Photons are an appealing infor- 
mation carrier for their inherently low- noise, high-speed 
transmission and the fact that entangling interactions 
between photons can be achieved using only linear op- 
tical circuits [TJ[T3-15 or mediated by atom-like systems 
[T6] [TT] . A photonics approach to these technologies re- 
quires complex, multi-port quantum circuits — essentially 
multi-path, multi-photon interferometers — that exhibit 
high fidelity quantum interference. Circuits fabricated 
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FIG. 1: Multimode interference devices, (a), Schematic rep- 
resentation of a 4x4 MMI integrated chip, (b) Simulation 
of classical light propagation in the device shown schemati- 
cally in (a). Light is launched into input waveguide 2 and 
multimode interference in the central region results in equal 
intensity in each of the four output waveguides, via self imag- 
ing. Analogous behaviour is observed for injection of light in 
each of the other input waveguides. 



from 2x2 directional couplers have demonstrated high 
performance [4-7, 18 , however, construction of more so- 
phisticated multi-port circuits would require their de- 
composition into a very large number of 2 x 2 direc- 
tional couplers (or X-couplers [Hj). For example, an ar- 
bitrary NxN mode unitary [20] would require a sequence 
of 0(n 2 ) individual 2x2 directional couplers. 

MMI devices are based on the self imaging principle, 
by which an input field profile is reproduced in single or 
multiple images at periodic intervals along the propaga- 
tion direction of a multi-mode waveguide [211 G2] • The 
effect is based on the propagation properties of a guide 
with a large number of lateral modes. MMI devices allow 
the design of NxM splitters with superior performances, 
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FIG. 2: Experimental setup for two photon quantum interference measurements in MMI devices. For the 2x2 MMI splitter 
filters A were 2 nm FWHM; a filter B of 0.5 nm FWFM was added to increase the coherence length. For the 4x4 MMI splitter 
two 0.5 nm interference filters were used as filter A and no filter B was present. 



excellent tolerance to polarization and wavelength varia- 
tions and relaxed fabrication requirements compared to 
the other main beam splitting technology, the directional 
couplers. Consequently MMI couplers have found appli- 
cations in a broad range of photonic systems [23] includ- 
ing phase diversity networks, light switching and modula- 
tors, in laser architectures and for optical sensing applica- 
tions. In the context of photonic quantum circuits they 
promise to dramatically reduce the complexity of such 
circuits, including for example those required to gener- 
ate maximally entangled path or 'NOON' states [24], W 
states [25] and the implementation of NxN unitaries [20 . 

In contrast to directional couplers, the self-imaging 
effect in MMIs allows the flexibility to directly realize 
symmetric NxN multi-port devices with several input 
and output ports. Multi-port circuits are particularly 
promising for quantum optics and information purposes 
and fundamental experiments have been conducted to 
study the behaviour of non-classical interference of sin- 
gle photons in bulk optics [26] and fibre [27] circuits. 
However, their performance is limited by stability and 
control of the splitting ratios. The implementation of 
multi-port splitters in MMI devices should allow higher 
performances due to the monolithic and scalable archi- 
tecture, with the caveat that quantum interference has 
not been observed in such devices: multimode proper- 
ties (see Fig. [lib)) might at first seem undesirable in a 
quantum circuit since they typically perturb or destroy 
quantum interference. 

MMI devices, including 2x2 and 4x4 couplers, were 
designed and simulated using a commercial beam prop- 
agation package (Fig. 1(b)). They were fabricated on 
a 4" silicon wafer, onto which a 16 /im lower cladding 
layer of thermally grown undoped silica was deposited, 
followed by a 3.5 /im core layer of silica doped with ger- 
manium and boron oxides deposited by flame hydrolysis. 
This core layer was patterned into 3.5 /im wide single- 
and 15 and 29 /im wide multi-mode waveguides via stan- 



dard optical lithographic techniques and then overgrown 
with a 16 /im upper cladding of phosphorus and boron 
doped silica with a refractive index matched to that of 
the lower cladding; simulations indicated single mode op- 
eration at 780 nm. The devices are composed of N single 
mode waveguides that serve as input and output for the 
multi-mode section and terminate with a separation of 
250 /im at the edges of the device to allow input and 
output coupling with a polarization maintaining fibre ar- 
ray. The 2x2 MMIs are composed of two input and two 
output waveguides that have a separation of 11 /im at 
the interface of the single- and multi-mode region. The 
multimode region measures 1090 x 15 /im. In the case of 
the 4x4 MMIs, the waveguides are separated by 8 /im at 
the interface to the multi-mode section, which measures 
1770 x 29 /im. 

Ideally, a balanced 2x2 MMI splitter should perform 
the same operation as a 2 x 2 directional coupler with a 
unitary matrix that describes the evolution from input 
to output [5]: 
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which equally superposes the two modes (equivalent to 
the Hadamard operation with a phase shift [35]). How- 
ever, it is not clear that the multi-mode nature of MMI 
devices will allow quantum operations, in particular 
quantum interference. 

Quantum interference with two photons is a defining 
distinction between classical and quantum states of light 
and is the key phenomenon that drives photonic quantum 
technologies. Quantum interference occurs when differ- 
ent quantum mechanical outcomes are indistinguishable. 
In the case of two photons entering the two input ports 
of a symmetric 2x2 unitary beamsplitter (one photon 
per input port) the outcomes of "both photons reflected" 
and "both photons transmitted" are indistinguishable. 
In this case the interference is destructive so that the 



photons never leave in the two separate outputs, but a 
superposition of two photons in each output. This be- 
haviour is in stark contrast to the case of two classical 
particles which would have a probability of 1/2 to leave 
in separate outputs. When the relative arrival time of 
the two photons is scanned, a characteristic "HOM" dip 
is observed [28] because the classical probability of 1/2 
holds for finite delay and the quantum probability of zero 
holds for zero delay. The width of this HOM dip is re- 
lated to the coherence time of the photons. The visibility 
V G [0, 1] of the dip (how close it gets to zero) is a mea- 
sure of the degree of quantum interference. Any informa- 
tion that distinguishes the two probability amplitudes — 
eg. the photons have different polarizations, frequency, 
bandwidth etc. — reduces V < 1. A key advantage of di- 
rectional coupler devices is that single mode waveguides 
mean that no spatial mode distinguishability is possible. 
In contrast, MMI devices are by definition highly multi- 
mode in the interaction region. 

We measured quantum interference in MMI devices 
using single photon pairs produced in the sponta- 
neous parametric down-conversion (SPDC) source shown 
schematically in Fig. [2] a type-I BiBO crystal pumped 
with a 60 mW CW laser diode at 402 nm, producing 804 
nm pairs of photons. These photons were collected into 
single mode polarization maintaining (PM) fibres after 
passing through 2 nm filters. The source was constructed 
in such a way that quantum interference with V ~ 98.5% 
was routinely observed, confirmed with a directional cou- 
pler that has previously exhibited V = 1 (Ref. [7j) [36] . 

We observed the HOM dip shown in Fig. |3(a)| in a 2 x 2 
MMI coupler. These data provide conclusive evidence 
that quantum interference does indeed occur in a MMI 
device (the linear slope in these data is due to decou- 
pling of the input fibre as the timing delay is changed) . 
However, the measured V = 90.4 ± 0.4% is significantly 
lower than the V « 98.5% obtainable from the SPDC 
source. The reason is that the propagation in the multi- 
mode section of the MMI introduces some distinguisha- 
bility between the photons. We experimentally ruled out 
spatial, spectral and polarization mismatch of the pho- 
tons, implicating the temporal degree of freedom. The 
different modes in the multimode section of the device 
see different effective refractive indices, which introduces 
a jitter in the time of flight of the photons from the in- 
put to the output waveguides, providing "which path" 
distinguishing information, and thereby reducing V. 

To confirm that this temporal jitter effect is the origin 
of the reduced visibility, we inserted a narrower 0.5 nm 
filter in one of the output modes between the device and 
the detector (as indicated in Fig. J2|, i.e. not affecting 
the properties of the photon source, but simply increas- 
ing the coherence length of the photons. The additional 
filter acts as a quantum eraser [29 , that erases the timing 
information by increasing the coherence time of the pho- 
tons. Under these experimental conditions, we observed 
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FIG. 3: Quantum interference in a 2 x 2 MMI coupler, (a) 
The measured HOM dip for 2 nm filters, corresponding to a 
dip FWHM of 239 fim. (b) The measured HOM dip for the 
same device and source, but with an additional 0.5 nm filter 
inserted into one output, resulting in a dip FWHM of 296 \ivn. 
The blue data show the measured rate of accidental counts. 



the HOM dip plotted in Fig. |3(b)| in the same 2x2 MMI 
device. In this case V = 95.6 ± 0.9%, confirming that 
timing jitter limits the visibility for the data shown in 



Fig. 3(b) (the larger error bar is due to the lower count 
rate with the narrower filter). These data confirm that 
quantum interference occurs in MMI devices, and that 
the coherence length of the photons must be sufficiently 
long compared to the timing jitter that is introduced as 
a result of the different refractive indices of the MMI 
modes. 

Interestingly, the description of multi-port splitters 
grows in complexity with N. For N < 3 all symmetric N 
x N splitters can be described by one equivalence class, 
since the requirement on the conservation of energy de- 
fines the matrix to within external phases on the input 
and outputs. However, when N > 4 there exists an in- 
finite number of distinct equivalence classes [30] and in- 
ternal free phases are independent of the conservation of 
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FIG. 4: Quantum interference in a 4 x 4 MMI coupler. Coincidence counts of two photons at the output ports of a 4 x 4 MMI 
device as the arrival time of the photons is varied. The different graphs represent the six possible input states in the splitter: 



(a) |ll) ia , (b) |11} 13 , (c) |11) 14 , 
interference filters used. 



(d) |11) 23 , (e) |11) 24 , (f) |11) 34 . The FWHM of - 800 \im is as expected for the 0.5 nm 



energy [24] . The transition matrix that describes an ideal 
symmetric 4x4 MMI splitter is: 
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where is the free internal phase. In general, two dif- 
ferent physical implementations would correspond to a 
different equivalence class, and a different value of 0. 

In the case of a MMI splitter the value of the internal 
phase is, in principle, dictated by the self-imaging con- 
dition [23] . However, the presence of fabrication imper- 
fections in the device would drive the multi-mode section 
away from exact self-imaging, and the relation between 
the optical phases would deviate from the expected value. 
Moreover, the presence of unavoidable losses in the MMI 
coupler corresponds to the presence of additional optical 
modes in the transition matrix of the splitter, thus mak- 
ing the reduced 4x4 transition matrix M more compli- 
cated than Eq. [2] In principle, M could be reconstructed 
using a number of phase sensitive measurements. Such 
measurements are, however, difficult in practice, due to 
the need to maintain sub-wavelength stability in inter- 
ferometers consisting of waveguide and fibre and/or free 
space paths. We have developed a technique to overcome 
this using only intensity measurements and two photon 



quantum measurements, but no phase sensitive measure- 
ments, as described below. 

In contrast to the 2x2 MMI device in which quan- 
tum interference is destructive, the 4x4 MMI device 
interference between indistinguishable outcomes can be 
constructive for some of the 36 possible input and out- 
put combinations [2 6 . We characterised the quantum 
operation of a 4 x 4 MMI splitter by inputing the state 
1 11) . . — i.e a single photon in each input waveguide i and 
j. We considered all six combinations of two photons in 
four inputs, where i ^ j (we did not consider the case of 
both photons in the same input since this does not give 
rise to quantum interference — i.e. such measurements 
provide no more information than bright intensity mea- 
surements do). As in the case of the HOM dip, quantum 
interference is revealed in the correlations in the output 
probability distribution. The probability to detect one 
photon in each output k and /, when two indistinguish- 
able photons are injected into inputs i and j is given by: 
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where Sij is Kronecker's delta and M^ is the element of 
the transition matrix. In the case of two distinguishable 
photons (equivalent to the classical analogue) the prob- 
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FIG. 5: Quantum interference matrix of a 4 x 4 MMI coupler. Measured (a) and reconstructed (b) visibility matrices Vijki of 
the non-classical interference between two photons injected in input waveguides i and j and detected in output waveguides k 
and / of a 4x4 MMI splitter. Positive visibilities correspond to a HOM-like dip, negative to peaks. 



ability is given by 
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and there is no interference between the two terms. Mea- 
surement of the generalized non-classical interference be- 
tween two photons enables the reconstruction of the ma- 
trix M via measurement of the detection probabilities 
Q\ l - and CtK as described below. 



Figure [4] shows the rate of detecting two photons at 
the six possible output pairs of waveguides, for each of 
the six different input states of a 4 x 4 MMI splitter, as a 
function of the relative arrival time of the photons. For 
some input-output combinations, these data exhibit clear 
interference dips, analogous to the HOM dip. However, 
for other combinations there are peaks, and essentially 
straight curves. These behaviours are the result of the 
phase values of the matrix M and are conveniently sum- 
marised by plotting the visibility of the non-classical peak 
or dip, given by: 
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where positive values indicate a dip and negative values 
a peak. Fig. p[a) shows the 6x6 matrix of measured 
visibilities V m obtained from the data of Fig. [4j 

We have developed a technique that uses only the val- 
ues Vijki and the classical intensity ratios |M^| to recon- 
struct the (reduced) transition matrix that describes the 
MMI device, assuming linearity of the device. To do this, 
we numerically search for a matrix M r that minimizes 
the RMS distance between the experimentally measured 
V m and the reconstructed V r that corresponds to M r . 
The classical intensity ratios \M^\ were measured using 
a CCD camera at the output of the integrated chip. The 
numerical optimization produces [37 
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The quality of the reconstructed transition matrix M r 
is confirmed by comparing the measured V m and calcu- 
lated V r matrices in Fig. [5] Additional optimization 
routines were performed with different initialization pa- 
rameters, to confirm that the measured transmissivities 
do indeed provide a good starting condition for the ma- 
trix M r . It is interesting to note that this method to 
reconstruct the transition matrix can be used for any 
unknown linear optical element (where no decoherence 
or measurement is present), even in the case of losses. It 
is possible to calculate the reduced n x m transition ma- 
trix even in the case of large NxN (N> n, m) networks 
in which only n input and m output modes are accessible 
for measurements. 

Interestingly, the computation of the photon distribu- 
tion at the output of a multi-port circuit with many pho- 
tons is a hard problem to solve from a computational 
perspective. The coincidence probability distribution is 
related to the permanent of the matrix that describes the 
multiport device [31 . Since the computation of the per- 
manent of a matrix is a computationally difficult prob- 
lem, the measurement of the coincidence probability in 
a device with many photons and many ports could rep- 
resents a feasible implementation of quantum simulation 
that uses the bosonic nature of the photons to realize a 
hard computation. In general, it is possible to map the 
calculation of the permanent to the detection of mul- 
tiphoton coincidences in the appropriate splitter. Al- 
though the form of the transition matrix that charac- 



terizes the MMI device is in principle fixed by the self- 
imaging condition, and no reconfigur ability is possible, 
the data presented here on a 4 x 4 MMI splitter is the 
first small-scale example of such a quantum computation 
[32]. 

Multi-photon inputs to multi-port devices are not only 
an essential ingredient of future photonic quantum tech- 
nologies, but enable the study of a rich variety of quan- 
tum interference phenomena to be studied [33 . The in- 
creasing need for more complex photonic networks will 
present the problem of the characterization of such inte- 
grated circuits and their imperfections. As demonstrated 
here, it is possible to take advantage of the properties of 
quantum states of light to reconstruct the behaviour of 
a photonic network, without the need of complex phase 
stable measurements. In the case of quantum networks, 
MMI devices promise to dramatically simplify and minia- 
turize the implementation of quantum circuits, thanks to 
the possibility of performing complex multi-mode evolu- 
tion of many photons in a single compact device. 
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